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KATIOKAL ADVISOEY COMMITTEE FOR AEEOMUTICS 
TECHNICAL NOTE 2506 

AN ANALYTIC DETERMINATION OF THE FLOW BEHIND 
A SYMMETRICAL CURVED SHOCK 
IN A UNIFOEM STREAM 
By C. C. Lin and S. F. Shen 

SUMMARY 


The method of power-series expansion in solving the local flow 
pattern behind a detached shock was proposed by Lin and Rubinov. In 
this report the limitations of the method are discussed in a more careful 
manner, and the practical procedure for approximating the power series 
with a 2nth-degree polynomial by cutting off the remaining terms is 
investigated. i 

It is pointed out that if the power-series expansion is to hold near 
the nose, the body shape must be analytic up to and including the sonic 
point. For a 2nth-degree polynomial approximation, n parameters deter- 
mining the shock shape, as well as the detached shock distance itself, 
are found to be expressible in terms of only n - 1 parameters deter- 
mining the body shape. The formulas and steps for a sixth-degree poly- 
nomial approximation are explicitly given in an appendix. 

The particular exanple of fr'ee-stream Mach number 1.7 in the axially 
symmetrical case has been worked out with a foTurth- degree polynomial 
approximation. When the detached distance is used as the length scale, 
the flow along the axis is found to be independent of the body shape in 
this approximation. A conparison is made with data obtained from an 
interferometric study of the flow over a sphere. The density variation 
along the axis agrees very well. The, detached distance as solved from 
the fourth-degree approximation, however, is correct only in the order 
of magnitude. 

The universal density* variation along the axis, as obtained by the 
fourth-degree approximation, is considered likely to be a good approxi- 
mation for a.1 1 bodies having a slowly varying curvature up to the sonic 
point at free -stream Mach nxxrabers larger than or not much less than, say, 
the value 1.7 considered above. 
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INTRODUCTION 


It is well-koown that a detached ciirved shock is formed in front 
of a hody vhen a supersonic stream flows past it, provided the hody has 
a hlunt nose or a sharp nose with a quite large nose angle. A general 
treatment of such flow problems is difficult. It is therefore thought 
useful to try to obtain some results of limited applicability by straight- 
forward methods of analysis. One such method is the use of power series 
proposed by one of the authors in reference 1 and partially carried out 
by Dxigundji in reference 2 . The present report is concerned with a 
closer Investigation of the validity of this method and includes a more 
coD5)lc be formulation, as well as the computation procedure, of the series 
expansion. General, formulas were derived and some numerical calcvilations 
were carried out. 

The case with free- stream Mach number I.7 has been confuted in the 
axlsymmetrical case using a fourth-degree polynomial to represent the 
flow. As stated in reference 1 , the flow thus determined depends only 
on the Mach number when the detached distance is adopted as the length 
scale. The density vaxiation along the axis shows remarkable agreement 
with the experimental one for a sphere (reference 3 )^ in spite of some 
indications of slow convergence of the series at the body nose. The 
effect of this poor convergence is manifest when an estimation of the 
detached distance itself is atteupted. The discrepancy is large, although 
the order of magnitude is correct. This fact leads to the suspicion that 
perhaps the variation of density distribution does not -depend very much 
on the exact shock shape and is therefore less sensitive to the inaccu- 
racy of the method, yet the detached distance itself is very delicate. 

It is clear that more terms of the series are needed for a better agree- 
ment of the detached distance. General formulas and the procedure of 
calculation are given in this report for such purposes. 

Because of the rather involved expressions for the series expansion 
and also the slow convergence occurring in the computed example, it is 
thou^t that the procedure here would be useful generally in the cases 
involving a body of approximately constant curvature in the neighbor- 
hood of the nose placed in a stream of Mach number not close to unity. 

One might notice that Busemann in reference 4 has stressed instead the 
iniportance of the "shoulder point," which presumably lies near the sonic 
point, in the determination of the detached shock in front of a body. 

His discussions, however, seem to deal mostly id.th thin bodies at rather 
low Mach numbers. A thin body with a blunt nose would have invariably 
a rather rapidly changing curvature near the nose. This configuration 
is then not directly amenable to the treatment by power series, (cf. 
discussions in the section "Discussion of Method.") 
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SYMBOLS 

stream function 

angle between incoming stream and shock wave 
ratio of specific heats of gas 

density of gas; also density ratio across shock wave 
pressure ratio across shock wave 

distance between shock wave and body along axis of 
symmetry 

dependence of entropy on streamline 
Mach number of incoming stream 
Increment of entropy across shock 

shock shape parameters, defined by equations ( 8 ) and 
(13a) for axi symmetrical and two-dimensional cases, 
respectively 

body shape parameters, defined by equation ( 2 l) 

functions in expansion for p, defined by equations (15) 

functions in expansion for defined by equations (l 5 ) 

functions in another expansion for ilf, defined by ■ 
equations ( 24 ) 

coefficients of Taylor series expansion of p and X, 
respectively, defined by equation (26) 

quantities evaluated at stagnation point 



k 


MCA TW 2506 


STATEMENT OF PROBLEM AND METHOD OF SOLUTION 


Consider a symmetrical body with a blunt nose placed symmetrically 
in a uniform supersonic stream. It is plausible from physical grounds 
to assume tbat^ for an analytic body shape, the detached shock vould 
likevise be analytic. 1 The Investigation will first be limited to such 
cases. This restriction evidently excludes some lniportant problems, such 
as the detached shock formed in front of a cone or wedge, which must be 
treated separately. 

To solve the problem, one may proceed in two steps. First, thi 
shock cur\'‘e may be regarded as given and tlie flow field regarded as 
determined by the initial conditions on the shock. If the shape of the 
shock is "reasonable," a stagnation point may be expected on the axis 
of symmetry, arid a streamline may be found in the form of a body. The 
analysis may then be carried out by calculating the flow behind an 
analytic shock in the form of a power series. Strictly speaking, the 
analysis does not depend at a1 1 on the presence of the’ body and consists 
only of finding the relation between the shock shape and the flow behind 
it. 


The second step is the investigation of the dependence of shock 
shape on body shape. Since the flow is determined by the shock curve 
and must also have the body as a streamline, the relation can be obtained 
by identifying the power-series solution from the shock with another 
power series developed around the nose of the body, the identification 
being made in a region common to the regions of convergence of both 
series. In fact, detailed calcxaations show that, by using series up 
to powers of the 2nth degree, this method enables one to solve for n 
"shock parameters," as well as the detached distance, in terms of n - 1 
"body parameters." The parameters are in fact chosen to be the curva- 
ture and its successive derivatives evaluated at the initial point (nose) 
Although the method appears straightforward, careful investigations are 
needed to insure its validity. These will be carried out in the section 
"Discussion of Method." 

The general analysis i- riTl be ceirried out in a manner applicable to 
both the two-dimensional and the axially symmetrical cases. Emphasis 
v;ill then be put on the latter case in working out the explicit formiilas 
and procedure. 


lln th.e actual case, there is a separation of the boimdary layer 
often occurring in the supersonic region. The downstream shock will, 
then cease to be analytic. But it is clear that this i-d.ll produce no 
effect on the flov; regime upstream of the separation point with which 
the present resort is concerned. 
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1116 basic equations for flow behind a symmetrical shock produced 
in a uniform stream are, by introducing the conventional stream 
functlo*n ill, 


y26 2p2(x ^ y ) + 


7F(^)p 


7-1 


7-1 


= C 


( 1 ) 


5/1 5itf \ ^ 5/1 5t \ 
5x\py^ 5x/ 5y\py^ 5y/ 


'(^) ye 


( 2 ) 


where x is along the. axis of symmetry, y is perpendicular to x, 
p is the density of the gas, 7 , the ratio of specific heats, and C, . 
Bernouilli's constant (the. same in front of and behind the shock), lie 
parameter 6 is zero for the two-dimensional case but equals unity for 
the axially symmetrical case, and F(ilr) gives the dependence of entropy 
on the streamline.. To reduce equation (l) to dimensionless form, the 
free-stream values of velocity and density may be chosen to be the 
respective reference quantities. The len^h scale is arbitrary for the 
moment. Then it follows that 

F(t) = (7M<„2)-leA8 


C = ± + ± 

2 (7 - 1)m„2 


(3) 


where M<j 3 is the free-stream Mach number and As is the change of 
entropy across the shock (in multiples of the specific heat at constant 
volume) . In fact, if a is the angle which the shock makes with the 
free stream. 


„Ab ^ I 27 
\7 + 1 


sin^a - 


7 - l\/7 - 1 


7 + l/\7 + 1 (7 + 1)M^^ sin^ 


(4) 
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The pressure ratio tJS across the shock is 


tS = 


27 

7 + 1 



sin a - 


7-1 

7 + 1 


( 5 ) 


and the density ratio across the shock is the dimensionless density 


P 


\7 + 1 (7 + 1)H»^ sin^a/ 


It is easily seen that 


e^ = 35P-7 


( 6 ) 


( 7 ) 


Althou^ formulas (i<-) to (6) hold only along the shock, the dependence 

of (7Mi^) ^S3p”^ on ■>{(■, as given hy equation (7)^ is valid throughout 

the whole field. One can therefore evaluate F(t) from the shape of the 
shock curve. 

Shock conditions.- Suppose now the shock cvurve is given hy 


2 

z = y*^ 




(8) 


Then it can. be easily verified that 

(sin a)"^ = 1 + kzy/[z'(x^^ 

= 1 + - 16S23i"^z2 - 

(33331"^ - 832^3 i "^) z 3 + . . . 
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(sin o ,)2 “ ^ 

= 1 -,l^^-2z + i6(p2 + l)Pi"^z2 + 

8Pi-2|3P23l-3 + 8(P2 + l)2pi-^lz3 + . . . 

By straightforward calculations from equations ( 5 )^ ( 6 ), and (7), one 
then obtains 
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where the subscript n denotes qusintities immediately behind ti^e nose 
of the shock. 


%i = 


2.7 

7 + 1 



7 -1 _i _ 7-1 2 

7 + 1 Pn "7 + 1 '^ (y + i)t^2 


I’n 



27 
+ 1 



2L^V2L_^ + 

7 + l/\7 + 1 


2 y ' 

(7 + l)^1co^/ 


( 10 ) 


and the coefficients c^, C 2 , c^, C^} 0^ depend only on 

the Mach number of the undisturbed stream: 


Cl = -87 - (7 - l)Moo"^ 


C 2 = -i»C3_ 


C3 = -2Ci 


Cl = 8/[(7 - 1)M„2 + 2] 


= -4c 1 


C 3 - - 2 c^ 


( 11 ) 


Equation (9) holds only on the shock. Hence, z may be replaced by a 
function of x through equation (8) . For exan 5 )le, one obtains 


p = pj^jl - C]^3]^“^x + {c-^ + 3c232)Pi”^^ ■ ’1 


( 12 ) 


All the above calculations apply to the conditions at the shock 
both for the two-dimensional case and for the case of axial symmetry. 
The transformation of the third relation In equation (9) into one 
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inTOlvlng t differs in the two cates. For the two-dimensional case, 
on the shock, ’ 


t = y 
= ^/i■ 


»l/2 


(13a) 


and hence 


- F^|l + ^7c^ + ^ j^(rc2 + C2)32 + 


^2 ^^1^1 


- 1) - 2 I 


- “3i-V + - C3[^3^P3 -8(32 + 1 )^ + 

7 CiC 2(32 + 1) + yCl(— 2 ^ ^ 1 ^ + ^ 232 ^ + 7^3(33133 - 8p2^) + 

r(7 - l)c,c,p, . c,3L^-6*6 , . . 


throughout the, whole field of flow behind the shock. In the case of 
axial symmetry. 


1' = iy2 


1 

= 2 ^ 


= §(3ix + 32^^ + 33X^ + . . .) 


(13b) 
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on the shocks and hence 


F(i{f) = 


Fj^jl + 2 ^ 7 Cq^ + ^ *^2)^2 

l^]|pi"H^ + 8|c3[33ip3 - 8 (P2 + 1 )^ + 




7(7 - 1 ) 


^^1*^2(^2 ^ ^^1^ — 2 — ^ ^2^2^ '*' ^^3 ( 33 x ^3 “ ^^2^) '*’ 

r(r - i)c:,c2P2 * ~ — P ^ ' -' *■ ■) U^i^) 


throughout the whole field of flow behind the shock. Equations (12)^ 
( 13 ), and ( 1 ^), with parameters defined by equations (lO) and (ll), are 
the basic relations associated irf-th the shock curve (8) . 


Differential equations .- To solve differential equations (l) and (2) 
in power series, it is found convenient to develop i and p first in 
powers of y with functions of x as coefficients. From symmetry 
considerations, one has 


^ = y^‘^^[^l(x) + t 2 (x)y^ + . . .J 

p = Pq(x) + Px(x)y2 + P2(x)y^ + . . . 


(15) 


Introducing these expressions into equations (l) and (2) and comparing 
powers of y, one obtains two series of differential equations for the 
two sets of functions ''Jfn(x) and p^(x). The function F(i|f) is given 

by equation (l4). 

For the actvial calculations, the two-dimensional case and the 
axially symmetrical case are best to be treated separately.' Only the 
axisymmetrical case will be carried out in detail. Again introducing 

the variable z = y^, one transforms equations (l) and (2) into 


2\ z 



+ 


7-1 


F(i)f)p’''*'l = Cp2 


(16) 
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P>' 


7-1 


Ft(ijf) 


(17) 


Equation (15) may be written as 


■>lf = ik2_(x)z + ■<li'2(x)z^ + i}f2(x)z3 + . . . 
p = Pq(x)z + P2(x)z^ + P2 (x)z3 + . . . 


(18) 


Introducing equation (l 4 b) in the abbreviated form 


F(ilf) = Uq + ajTjf + + . 


( 19 ) 


and substituting equations (18) and (l9) into equations (16) and (17) 
one derives the two sets of differential equations: 


N) 


‘‘V + 7 ^ ' ^<^0 






= “KPoPi 


(Aj) 2 ti'i 2 ' + ^(^ 2 ^ + p, 


7+1 

0 1 


(7 + 1 )^ + 


7(7 + 1) /Pi 


'2 VpQ 


+ + 1 )^ + a 2^2 + 

Pq 


^2^1' 


I = 2 C (pj^ + 2PQP2) 
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(^ 3 ) (^ 2 ' ) ^ ' ^3 ' + 4 {^2^3 + Pq^+1. 




(7 + l)7p + 
PQ 


7(7 + 1)^ + 2( 7^ -1) 


Po" 


■a^ti 


(7 + 1)— + 2(7 

PO 2 p^2 


(^1^2 + «2^i^) (r + 1)^ + a^t3 + + a3^^3 


h ^(PqPs + PiPs) 


N) 


i,f^K - £iU .-^ „_r 


dxlp^ 


P0\ti Po/ 7-1 ^0 


W s 


fi 

L^O \«i' ■ Po 


V3% Pl2«i Pi^ 


PQ^ / Pt \ 

T^fa2% + ?'ai 


H ^ 

Pi '^ 2 ' Pi® P2\ 

' dx 

Po W' ' Po ♦!' " Po® ■ Po/_ 


12 


^-4 — - 3 ^i + 2I— 23 

oL"i V'Poi^ 


•y 

■~^ 282+2 + 383+5^2 + 2782+, ^ + a, fr ^ + 2l(y - Ij 

L ^ ^ ^ Po H ‘’0 2 p^2 
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Equations and are obtained by equating coefficients of zP 

in equations (16) and ( 17 )^ respectively. The other equations arise 
from coefficients of corresponding powers. 

One notes tliat from equation (Aq), can be expressed in terms 

of pQ. Then from equations ^Aq^) and (^ q)j ^2 Pi expressed 

in terms of Pq, Pq', and Pq”; from equations ^A^) and 
P2 can be expressed in terms of Pq, Pq', Pq", Po"'-> ' PQ^^-' ^ 

Determination of solution from the shock curve .- It will now be 
explicitly shown how the solution can be calculated as a power series 
when an analytic shock curve is given. For definiteness, the discussion 
will be limited to expansion up to the’ following degrees: 


Pq(x), ■' 1 ''q^(x) up to the fourth degree 
Pp(x), ■'^■2(x) up to the second degree 
P2(^)j ^ constant terms only 

Referring to equation (I8), one sees that the physical variables p, 
u, and V are expanded to a foTirth-degree polynomial in x and y. 

In general, there is a decrease of two degrees when one proceeds 
from ^p^, \+l) (^n+1^ ^1+2) * ' directly related to the above 

discussion about expressing them in terms of Pq(x) and its derivatives. 

Again, in doing so, the expansion of p, u, and v as power series of 
X and y Includes all the terms up to a certain degree. 

Returning to the case of the fourth degree, one has to determine 
five coefficients for each of Pq(x), ilr^Cx), three coefficients for 

each of pj^(x), ilfgCx), and the two quantities p2(0) and 

However, since the quantities ilr^, p^, p^, and expressible 

in terms of Pq(x) and its derivatives (up to the proper orders), there 

are only five coefficients to determine. This is done by' imposing 
conditions ( 12 ) and (l 3 b) for p and f on the shock. These power 
series are compared with another form of the same development obtained 
by Introducing equation (8) into equation (18); namely. 
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P = Pq( 0) + 


p2(0)Pi^ 


Pq'(O) P]_P 2^(0) 


+ . . . 


X + 


I Pq"( 0) + Pi’(0)p3 ^ + Pi( 0)32 + 


ilr - I z = 




x2 + 


>|fl*(0)32 + I 


i2'(0)Pl^ + 2^2(0)3132 + 


x-* + 


( 20 ) 


Exactly five conditions are obtained involving only the quantities 
desired. The only geometrical parameters of the shock entering the 
problem here are 3i and 32 which appear in the above equations and 

also enter through Sq, ai^ and a2 appearing in the differential 

equations. Thus, the calculations of the solutions up to the fourth 
degree Involve only two geometrical parameters from the shock. Thi s has 
already been discussed by T i l .n and Rubinov (reference 1, p. 122), but 
here the procedirre is explicitly outlined. 

As a general rule, for a 2nth-degree expansion of the flow, the 
quantities Involved wi l l be up to '^^^^(O) and pj^(O). There are con- 
sequently n conditions from ilr and n + 1 conditions from p by 
inposing equation (20). Thus the 2n + 1 coefficients in the expansion 
can be \iniquely detennined . It is also noted that there are n shock 
parameters 3i, 32^ • • 3^ entering the problem. 

Relation between body shape and shock curve . - One starts next to 
expand the solution about the stagnation point of the body. Except for 
a shift of the origin, the previous equations (16) to (19)^ together 
with the consequences in the form of the equations (a) and (B), evidently 
hold just as well in the present expansion. The observation that 

and can be expressed in terms of Pq, Pq', Pq", . . ., Pq^^^ 

remains also true. Suppose now the body shape is given by 


z 



2 

6q^x + ^ 2 ^ + . . . 


( 21 ) 
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Then the conditlone that furnish the relations for the determination of 
the coefficients become 


Pq(0) 


/ 27 c 

^7 - 1 


0 


(22) 


and (cf. equations (20)), 


0 = 


ti'( 0 )ei + ^ 2 ^ 0 ) 


x2 


il'(0)€2 + I + 


^ 2 '(o)^y + 2 ^ 2 ( 0 ) €^ 6 ^ 



(23) 


One might reason that In comparison with equations (20), the con- 
ditions from density variations along the body are misBlng, and there- 
fore there are n + 1 fewer conditions. Even with the one for Pq(0) 

in equations ( 22 ), it seems that for 2 n + 1 coefficients in the exj)an- 
sion only n + 1 conditions are available, leaving a total of n arbi- 
trary parameters. Although the resxxlt turns out to be true, the argument 
is not so simple. An important point is that instead of, n body param- 
eters according to the above argument, only n - 1 parameters defining 
the body here enter the problem. This may be verified as follows. 
Rewriting 


ilr = yl+^X(x,y2) 

= + . . .) {2k) 


with X( 0 , 0 ) =0 because the origin is chosen at the stagnation point, 
one may transform equation (16) into 


. 7-1 ^ 


- ^ + 2(1 + €)yXyX + X^ 


7-1 


F(t) 


7 


(25) 
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Each "term In "the bracket of the right-hand side is seen to have a double 
zero at the origin. After successive differentiation^ the highest deriv- 
ative of X on the right-hand side appears to be at least one order 
lower than the highest order of p on -the left-hand side. In other 
words, if one writes formally 


P 


i J 


PijX-5yi 


(26) 


and similar expressions for the other variables, there follows 


p = cfx ,X , . . 

^ij i-2,j+l" i,J-l' j 


(27) 


where the omitted terms are lower-order derivatives of X and p. Since 
p is even in y, the index i is always even and terms like 

do not exist. Needless to say, the subscripts can never become negative. 
Men a representation, using a 2nth-degree polynomial is made, all the 
p . , ' s up to i + j = 2n are required, involving the following deriv- 

1 J 

atives of X: 




Xqi Xq2 . . . Xq 

,2n-l 

^1 • • -^,2n-3 


(28) 


^-4,0 ^^-4,1 • • • 

\ 

^-2,0 ^^-2,1 


a total of n^ + n - 1 qiiantlties. Among them eqimtion (l?) furnishes 
some relations. By means of transformation (24), it becomes 


’Sex + V - Mpx’Sc + PyV) 


bi _ (1 


(1 - e2). 


X 


p7+lF.(,jf) ^ ^ 


,1+e _ i)yl-e 


+ 


(29) 
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Then X. , is obtainable by applying the operator — ^ — r—^ to equa- 
ij 

tion (29). All the terms except the first two can readily be seen to 
consist of lower-order derivatives than xy, the' derivatives of p 

being representable by lower-order derivatives as expressed in equa- 
tion (27)- Hence 

Xij = H(Xi, ■2 j+2^ lower-order derivatives^ 

As a result, all the X^^j's are expressible in terms of the derivatives 
with respect to x only. A total of 2n - 1 parameters ^e.g., Xq]_^ 

Xq 2> • • ^0 2n-l) sufficient for the determination of a].l the 

needed ® s. 2nth-degree polynomial representation for p. 

The general expression of condition (23) is now 

-"*^2 51 + \{YL + . . . = 0 (30) 

, n=l \n=l / 

The coefficient associated with x™ will yield an equation involving 
X2 jjj^ 0 highest-order term of X and 6]jj as the highest-order 

term of body shape. Equation (28) indicates that f or a ' 2nth-degree 
expansion the highest-order term of the form Xp^ ^ needed is ^2n-2 O' 

Hence, equation (30) will furnish useful conditions by equating the 

coefficients to zero for terms up to ’ x°-~^, that is, a total of. n - 1 
conditions. There are thus n arbitrary parameters (as suspected) left 
in the determination of the expansion. But since’ only terms up to 

xU-1 used in equation (30), there are involved but n - 1 parameters 

defining the body shape. ' - 

It remains now to show how the relation between the shapes of thfe 
shock curve and of the body may be found. As stated above, the expan- 
sions from the shock and from the body are matched at a point ly ing within 
the regions of convergence of both. In particular, to avoid the, question 
of s in gularity which might occur inside the body contour, - the stagnation 
point itself is arbitrarily chosen as the "matching point." Along- the 
x-axis, the 2nth-degree expeinsion furnishes ’ 2n + 1 conditions on ' 
equating the successive derivatives at the. matching point. , The .variables 
Involved are: n parameters defining the shock curve 

Pn), - parameters in the expansion from the body (-e,g., p, ■ 

. . ., n - 1 peirameters defining the shape of the body ^ e£, 
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. . 1) distance 5 between the body and the detached 

shock wave along the x-axis. Hence a total of 2n + 1 unknowns may be 
solved in terms of the n - 1 parameters defining the body shape. It 
may be noted that the solution obviously would be modified when different 
values of n are taken. However, if the series expansions do converge, 
the result should tend to a limit for increasing n. 

If one does not care for a complete expansion in the form of eq\ia- 
tion (26) but rather is interested in the distribution along, say, the 
x-axis, the process is sometimes much simpler. By writing 

00 = PlO + PoiKf) + * ■ ■ ■* Po, 2 n^^“(f)“ (31) 


the coefficients Pq involve but a total of n parameters 

pg, • • Pq^)j with 5 as the length scale. One could pick out from 

the set of 2n + 1 simultaneous equations a subset whijCh is sufficient 
for the solution of P^j Pp^ • • •> Pn* general the subset would 

contain parameters occurring in the expansion from the body nose. In 
the particular exanple of a foirrth-degree polynomial, the two necessary 
equations for Pj^ and Pg are readily provided by equation (22) or its 

equivalent and the expansion from the body is not needed. Evidently the 
body shape can not affect the fourth-degree expansion so determined. 

The simplicity of this result is certainly too attractive to be ignored, 
in spite of its possibly poor accuracy in extreme cases. The result 
ttirns out that a good agreement in density distribution is obtained with 
the experiment over a sphere at Mo<> = 1 «T* Indication will be made in 
the section "Discussion of Method" as to the possible range over which 
a similar agreement might be expected to hold. When one tries to obtain 
the value of 6 by imposing other conditions in the matching process, 
however, the result for the fourth-degree expansion shows a large dis- 
crepancy in comparison with the same experiment quoted above. This fact 
is not surprising since the power series is cut off after only four terms 
and the error would be undoubtedly laxge when higher derivatives are 
taken. 


Explicit formulas and procedure .- In the actual carrying out of the 
solution, equations (A) and (B) are to be expanded further into ascending 
powers of x. Equations expressing the relation among the quantities 
Plj ^ij obtained by equating the coefficients of each power 


of X. These formulas and also the procedure to solve for the desired 
coefficients in the case of a sixth-degree polynomial representation of 
the flow along the x-axls of an axially symmetrical body are given in 
the appendix . 
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DISCUSSION OF METHOD 


The method described above^ though seemingly straightforward, cannot 
be relied on without a more careful examination as to its limitations. 

Two points must be considered. In the first place, the subsonic region 
just behind a detached s-hock is expected to depend on the entire body, 
shape. This seems to make it difficult to treat the problem by power 
series, which depend on local properties. The second is the question of 
convergence of the series. 

The first point is really the following question: To what extent 

is the flow field determined by the local properties of the body at the 
nose? Theoretically, if the body shape is analytic, the local properties 
of the body do determine the entire body shape and thereby determine the 
entire- flow field including the shock. In fact, at some distance down- 
stream there is a sonic line (CD in fig. l), after which the flow becomes 
supersonic. For two-dimensional flows, Guderley has suggested (refer- 
ence 5) that, based on Trlcoml's study of a certain differential equation 
which is elliptic in part of the region and hyperbolic in the rest, there 
should be also a unique solution for the subsonic region behind the shock. 
He reasoned that the body shape should be given up to the point E, from 
which the Mach wave reaches the sonic point on the shock front. Any 
reasonable modification of the body do-vmstream from this point -will not 
influence the flow in the region ABCED. The same conclusion, when 
physically incerpreted, seems to be squalls'^ valid in three-dimensional 
flows. Thus, if a s\ifficient number of terms are retained in the power 
series for the body shape to represent it over the part BCE, slightly 
beyond the sonic point, all the important parameters of the problem are 
known. If the body has a fairly blunt nose, a few terms would be suffi- 
cient. The method of series expansion may then be expected to work. 

It is iii5)ortajit to notice that the above statement applies only to 
cases where the body contour is analytic over BCE. If the sonic point 
were brought in by the presence of a sharp turn in the boundary, a 
Prandtl-Meyer expansion exists, and the flow field in the region ABCD 
woxild be determined in an entirely different manner. 

In order to apply the above discussions to a given body, one must 
be able to estimate the location of the sonic point. In certain cases, 
this might be done, for instance, by the method sxiggested by Busemann 
for the determination, of the shoulder point (reference k) . One might 
even suggest that the body shape should be more closely approximated in 
that neighborhood at the expanse of it’s nose. The important point is 
that, since the entire contour BCE affects the shock, certain judgment 
is needed in picking out the most significant parameter. When there is 
more or less constant curvature starting from the nose, and the detached 
distance is kno-wn to be small in comparison -with the radius of cur-vature 
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at the nose^ vtndoubtedly the nose curvatiare plays a predominant role. 

Local properties of the shock, such as the curvature at its nose, would 
he largely determined by the nose curvature of the body. On the other 
hand, in cases like those discussed in reference L, one has a thin body 
v;ith a blunt nose and a quite large detached distance. The body curvature 
changes at a considerable rate in the neighborhood of the nose. This 
i-ate, then, might have an effect of the same order as or even overshad- 
owing that of the curvature itself. The flow near the shock, being far 
away from the nose, would depend on the over-all body shape BCE. The 
simple procedure of using only a few terms of the expansion depending 
only on the body nose curvature is obviously Insufficient for any 
accuracy at all. 

After liaving clarified the dependence of the flow on the body shape 
near the nose, one must take up the question of the convergence of the 
series. Because of the analytic nature of the flow in the subsonic 
region, it might be assiuned that the power-series expansion from the 
point A on the shock is convergent up to the body, including the stag- 
nation point B. The region of convergence of the series starting from 
B, on the other hand, depends on the nearest singularity lying within 
the body contour. In choosing the point common to the regions of con- 
vergence for matching the two series, one must try to be close to the 
point B. In the section "Statement of Problem and Method of Solution" 
the point has in fact been taken to be point B itself. Indeed, if 
the series from A and B are closely approximated by polynomials 
of the same degree, the matching of the two series gives completely 
identical resirlts, no matter where the conditions are applied. In fact, 
the present procedure of applying the method may also be regarded as one 
of polynomial approximation, satisfying a finite nuinber of conditions at 
the boundaries . Then Instead of talking about convergence, one can say 
that the inaccuracy is due to the faJLlure to satisfy all the conditions 
as required by the governing equations of motion. There are techniques 
for improving the (over -all) accuracy of approximations of this nature. 

It is likely that by utilizing some of those techniques (e.g., least- 
square error) the practical value of the present method could be greatly 
enhanced. 

After assuming the convergence of the series representation of the 
flow variables, one still has to consider the error caused by the retain- 
ment of only a finite n\imber of terms. It is lflg>ortant to point out 
that the solution for the parameters changes with the degree of the 
polynomial chosen, although theoretically the values would converge as 
n increases. The appropriate degree n’ for a satisfactory result, 
therefore, can not be stated a priori. One could merely argue that, for 
a body shape with slowly varying curvature up to and slightly beyond 
tlie sonic point, the shock-wave curvature presiimably would also vary 
rather slowly up to its sonic point and, consequently, a polynomial of 
a fairly low degree would be needed in such cases. The effect of the' 
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Mach number on the convergence of the series is more difficult to 
visualize. The question strictly can only be settled through actual 
calculation. 


To clarify this point further, the validity of the fourth-degree 
polynomial as the universal function for density variatlon-between the 
detached shock and the body nose will now be discussed. As pointed out 
above, one must look into the neglected terms in the series expansion. 
Consider, say, a sixth-degree representation. The coefficients thus 


solved will now contain, besides the Mach number, the product 

where is the body curvatiu*e and 6, the detached distance. VJhen 


terms, 

small. 


is small, the modification in the coefficients of the first five 
as compared with the fourth-degree representation, would also be 

By dimensional reasoning obviously depends only on the 

Mach number. So a range of Mach numbers prevails over which the fourth- 
degree representation is very close to the first five terms of the sixth- 
degree one. A thorough investigation should include the variation of 
the coefficient of the last two terms in the sixth-degree representation, 
which would lead to another restriction on the Mach number to justify 
their omission. A common range of Mach numbers would presumably become 
available for the fourth-degree representation to hold. These^ steps, 
though desirable, were not taken in this report. Without going into 

details, it seems physically likely that the sma lln ess of would 

ensure the fourth-degree representation to be a reasonable one. For, 

the smallness of means that, for a given body curvature 

8 must be small. A smaller 8 very likely tends to improve the practi- 
cal convergence of the power series. It should be noted that, in char- 
acterizing the body by a single parameter one implies that the body 

shape must not deviate appreciably from a parabolic one up to the sonic 
point. 


If the criterion of small K^®^8 turns out to be correct, evidently 
a larger Mach number is favorable for the approximation. Semleng)irically, 

the smallness of K^*^^8 may be tentatively measured against the value 
in the experiment quoted above, to be discussed in the following section, 
where the measured density distribution agrees well with the fourth-degree 
representation. 


NUMERICAL RESULTS OF FOURTH-DEGREE POLYNOMIAL APPROXIMATION 

The c§.se of free-stream Mach number I .7 has been computed for an 
axially symmetrical body by using a fourth-degree polynomial as the 
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approximation to the power-series solution. The results indicate that 
calcinations Involving higher-order terms are desirable for the deter- 
mination of the distance of detachment. The steps outlined in the 
appendix ha'^e been followed. Taking 7 = 1.405^ one may list the results 
as follows: 


Poi = 7.296P1-1 

= -2.38k&^-^ 

“10 = -12-833 i -2 

11^20 = 2.38Upj^-2 

ilf ^2 = -16.51Pi"2 
2q2 = - 28-60|3^-2 
Pll = -29.l8P2pi“3 - i42.4pi-3 
= 4.768P2p^-3 + 53.63P^-3 
^30 = -7.152P2pi"^ - 37.l2Pi-^ 
P20 = 58.6ip2pi"^ + i 84.9 Pi'^ 

= -lo.8op2pi‘3 - 73 .^Pi ‘3 
P03 = 33.o4P2Pi‘3 + 390 . 7 Pi '3 
1^22 = l4o.4P2Pi"^ + 596 . 2Pi"^ 
P 12 = -344.iP2Pi“^ - l772Pi"^ 

= -il9.7P2Pi'^ - 550 . oPi"^ 
Pq4 = - 350.3P23 i“^ - 68o6p^-^ 


There actually was another set of these functions because Pq^^ is 

solved from a quadratic equation. However^ the other expression gives 
a negative value of Pqi^ contrary to the conception of a compression 

along the axis and is therefore abandoned. 
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Conditions (22) may now be used to solve simultaneously for the 
parameters Pp/5 £ind ^2* eqiiations are: 


0.479 = Pq-^S + Po2 ^^ P03^^ Pq 4^^ (32a) 

0 = 0.500 + ^ ^^^ q 3 + (33^) 

After substitution eqiiations (32a) and (33a) become 

o.i.79 = 7.296 ^ - 28.6o(^)^ + * 

(33. ou - 350.3 ^){^)\ (32b) 

0.500 = 2.384 ^ + lfi-5l(|;;)^ -b 73.116(^)2 + * 

(10.80 + 119.7 (33b) 


By eliminating P2.> ^ fifth-degree algebraic equation in B/Pp is 

obtained. Again one has to choose the proper root which corresponds to 
the physical problem. . In the present case^ the choice is made easy by 
comparing with the experiment over a sphere as presented in reference 3- 
The proper root is found to be 


^ = 0.0644 
Pi 


P2 = 50.4 


( 34 ) 
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Consequently, the density variation along the axis 


as 


Po = 2.192 + o . k 699 § - 0.1186(1)^ + 0.5^93(1)^ - 


may be written down 

0 . 4214 ( 1 )^ ( 35 ) 


Up to this degree of representation, the body shape does not enter. 

Before making comparison with the experimental data in reference 3^ 
it may be noted that the solution for 5/Pi and ^2^ such as equa- 
tion (34), depends on the conditions chosen for their determination. To 
illustrate this point, one may recall that, instead of the condition 
= 0 at the stagnation point, an equivalent condition is dpQ/dx = 0 

at the stagnation point along the axis. In the simultaneous equation, 
equation (33a) is to be replaced by 

0 = PqiS + 2po25^ + 3 Pq38^ + 4pq 4&^ (38a) 


so that in parallel vrLth equation (33t) one now has 


0 


= 7.296 ^ - 57.2l(i^)^ + * 

(99.11 - 1401 


(36b) 


The solution of equations (32b) and (36b) gives 


= 0.0697 
Pi 

P2 = ^T-5 


( 37 ) 


as well as 

Pq = 2.192 + 0.5085 I - 0.1390(1)^ + 0.6633(1)^ - 0.5538(1)^ (38) 
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The difference between the two sets of results is not su:^rising, because 
only a finite number of terms have been tised in the series. The- terms , 
omitted have imequal effects on the approximations for Pq, dpQ/dx, and 

Equations (35) and (38) are .both plotted in figure 2 together with 
the experimental points from reference 3- The curve for equation (35) 
gives a slightly smaller density Increase than- equation (38) throughout 
the reinge. The agreement with experiment is good except for points very 
near to the body. There the experimenthL value is much higher than the 
theoretical one assuming isentropic compression. The discrepancy is too 
large to be attributed^ say^ to the neglect of the viscosity effect in 
the theory. On the other hand, one may suspect that the accuracy of the 
experiment could be poorer in the neighborhood of the stagnation point, 
where a very small region of high density occurs (see fig. 29^ refer- 
ence 3)* A slight misalinement of the apparatus with the instantaneous 
actqal flow direction, for Instance, may cause a quite noticeable change 
in the value evaluated from the photograph. 


The condition dpo/dx is, in fact, the one which results from the 
genereil scheme- of matching the expansion from the shock with the expan- 
sion from the nose. Such a scheme will now be worked out for the present 
example for the purpose of finding the, detached shock distance 6. ■ Let' 
p and if denote the variables in the expansion from the body nose. - • 
Again, by following the steps outlined in the appendix, one obtains 
readilyi 


Poo 


1 

y - 1 
^ ^0/ 





^20 

tl2 

P02 

P03 


1 * 2 




11 

-1 


Dl 

-^11^1 
%1«1~^ 

El 


- y+l 

2(7 - 1) ^EO 


11 


8 T 


Huei 


-1 


- 7+1 

2(7 - 1) 





- ' 8 T T 'V 

Pol. = ^11^13 - 


^^ll"l 


-1 


^1 - v+l 

2(7 - 1) *^00 


E2 _ 

” Dl ^ 


02 
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The matching then involves a £et of five slmulteoieous equations for the 
five unknowns P 2 -> ^U-> ^ terms of the body sliape 

parameter 


PQI^ ■" + 3P0363 + 4 Pqj^ 6^ 

2Pq2^^ + ^Pq 3^^ ■*■ 12pQij,6^ 
^PQ3^^ + 24pqJ^5^ 

Po 4 &^ 


"00 


Poi^ 


2po2B^ 




6 pQ 363 

Po 45 ^ 


(39) 


Obviously the first two are equations ( 32 a) and ^ 3 ^a) solved before for 
the parameters S/Pp anh &2* Note also that 1^23 appears only in 

pQlj^j therefore the first foiir of equations (39) are sufficient to deter 

mine P-j_, Pg^ 5 , and easily verify that 


^1 SpQgS^ 


2(r - l)Di Pog62 


(4o) 


The negative root is taken for 

because ~ ^ at the stagnation pointy which is negative in the 

physical problem. With solution ( 37 ) 7 "the detached shock distance turns 
out to be 


with 


" -(-X 


= 0.138 ( 4 l) 

^1 

In comparison with the data in reference 3 ^ the value of 6^ being 
equal to the diameter of the sphere, the approximate formula ( 4 l) gives 


6 = 1.1 mm 



KACA. TN 2506 


27 


while the experimental value is 1.9^ millimeters. The discrepancy is 
certainly large^ yet the order of magnitude is correct. Since only a 
polynomial of the fourth degree is used, giving a density distribution 
independent of body shape, this is perhaps the reasonable extent of 
agreement that can be expected. More careful calculations must Involve 
terms of higher orders . • 

One may also be interested to see how the shape of the shock curve 
is being approximated by the above solution. With equations (37)^ the 
shock shape from the present approximation is 

(if '= • m 


The curve (42) has been compared with the experimental one given in, 
reference 3* The agreement is reasonable for small values of y/S, but 
the deviation becomes appreciable for values of y/S of the order of 
unity, as may be expected from the slow convergence of the series (3^) 
and ( 38 ) . . 


Massachusetts Institute of Teclmology 
Cambridge, Mass., June 30, 1949 
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APPENDIX 

EXPLICIT EOEMULAS AND PROCEDURE EOR SIXTH-DEGREE 
POEXNOMIAL APPROXIMATION OF FLOW ALONG 
AXIS OF AXIALLX SYMMETRICAL BODY 


It is desired now to find an approximate solution in the form 


Po = Poo Poi^ + Po2^^ P03^^ Pol^^^ P05^^ ^ Po6^^ 


for the flow along the axis of ar. axial ly symmetrical body. As stated 
in the section "Statement of Problem and Method of Solution," 
equations (a) and (B) are to be ei^anded into ascending powers of x. 

The resulting equation from eqiiating the coefficients of a certain power 
of X of equation (Aq), say, is to be denoted by a second subscript. 

For instance, (Aqq) stands for the equation obtained by equating the 

coefficient of x^ in equation ^Aq) . In a similar manner, equation (B22) 

stands for the equation obtained by equating the coefficient of x2 in 
equation (B2)^ and. so forth. By so doing one gets the following: 


(^) 

An identity between = l /2 

Poo “ 

(^1) 



(^2) 

» 1»02 = -‘*(2+io'''i2 + 

- D2 Poi^ 

(^3) 

D1P03 ='-®{Vi 3 + ♦n% 2 ) - 

2D2P01P02 " 


27 

7-1 


„ f, 7 +I 0 7-2„ 3 

PQO Pq! 


O 
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(AqI^) . Vo4 “ "^(^10^14 ^11^13 \^) - {^W03 ■*■ PQ2^)®2 " 




2y 

7-1 




00 


(^ 5 ) %Po5 " (^ 10^15 ^ 12 ^ 13 ), -‘ 

(2poiPoI<. +• 2po2Po3)®2 ■ 


7-1 


^0 [^^3^^^^00^ ^(^01^2q 3 +• PoiPo2^) 
-7+1 7-3_ 3- - p 7+1 7-h 5 I 

^4 *^00 ^or ^02 ^5 ^00 *^01 J 


+ 


(<V)6) '>lPo6 = -*t(2<'io’l'l6'+ 2^11%, + 211-12+14 + +;^ 32 ) . 

(^'^0l‘’05 * ^Oh’‘o2 * '’ 03^)®2 ■ 

a^[c3)'+lpog+-2(3Poi2po4 -f Sp^^^p^jPcj + Pqs^) 

c/*\/-^(5Pi,iV) . Cg^-ip„/-5p„/1 


(%o) * 11 ^ + lfi* 10*20 + EiPio = - S.^00*^W 
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(^n) ^^(^10^21 ^20^11) + %Pll = 


2y ( *y *x^l 

■12'U —2'^10"'01 ' ^iri ^00 PqAo f^oo ’t'll 


- 2D,P,^P, 


(%2) V 12 = -®'''ia’''l3 - '*'''12^ - ^^{’I'l0'''22 + ♦li'''21 ■" ♦l2'l'2o) ‘ 

®2('^lo'’o2 * ’’nf’oi) ■ 7-^1 


27 

7 - 1 ^ 


’^i2‘’oo’'"^ + '=/*'''=oo V’^ia + 


(=l^*"P0O^P02 - =2^%0^- 


V^)’l'lO 


(^13) ® 1^13 “ “^ 11^14 " 12^12^13 " ^^(^10^23 + ^11^22 + 

^12^21 ^13^20) " ^2(f^loP03 ■*■ ^02*^11 ■*■ ^01^12) " 

1 ^oj|^3^ ^00^ ^(^PoiP02PlO 3 PiiPoi^) + 

"'l’"'''P00''P01^2 K’'"'''Poo’'Po2 ^ 


„ 7+1 7 7+1 7-1 7+1 7-2 

^1 ^00 *^03 ^2 ^00 ^01*^02 ■’■ ^3 P(X) PQl /^lO 
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(Ai 4) Vl4 = - ^13^ - ^Hho*2k * V23 + 


^12'*'22 ^13^21 ’''14^20) ^2(^01^13 ^02^12 


P03Pn 


11 '*' '^04^10) “ 7 - 1 ^0 ^3^ ^00^ (^^lo'^Ol'^03 ^ 


^“o2\o * ^0l“o^ll * 30o1^Pi2 I 




c/'^^P ’’•3 
4 ^00 


(^^^01 ^02^10 ^^01^^11 ) 


7+1 7-h 

^5 ^00 


C ijf + (c ^ 

1 ^00 *^01^3 ri 




^■*■^0 ’'p + C ’'+2 q 7- 

^00 ^02 ^2 *^00 


■VI 


t + 
12 


’'‘"^Pnn’'pnV+ ^o^'*'\r7~^Pn.Pr.^ + 


(^1’ *^00 *^03 "^ ^^2 »"00 ^0lP02 

^3^ ^^00^ ^^01^)% ■*■ [^l^'^^Poo’^Po4 + 

^2^ ^PqO^ ^(^^01^03 ^02^) ^3’^^^^00’'"^(^P01^P02) 


C,7+lp 7-3p 4 ]^ . 

^00 ^01 J \o 


(A 20 ) 2+3^iV + Vo»30 + ”l‘=20 * 

^ VOO^K’^Vo'k^g + PjQtjg) 


-Vo -Vio - 7 ^ 


27 „ 7+V 2 

7-1 ^2P00 UO 
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( A21 ) “ “^(^ 11^22 ■*■ ^ 2 ^ 21 ) " ^^^ 20^21 

^2 (*^ 20^01 ^ lo '^ n ) " 


_ pkilf ilr 

11 30 


7 ^ aoC 3 ’""^ Poq ^'^( 3 PoiPlO ^) - 
Pqo ^( ‘^ 10^11 ^lAo 

^ Ol '*^ 2 o ) ^ 2 ^ ^ 00 ^ ‘^ 10^01 ■*■ ^ 00 ^ ^ 2 ^ 


y~ Jl ag (2p00^’^^^10^11, ^1^ ^PQO^Poi^lO . ) 

1 

(^22) ®lP22 = -^^11^23 ■ ^^12^22 ■ ^’1^13’*^21 " 1^(2^20^22 ■*■ ^21^) " 

2^(^10%2 + ^U%1 + ^12%o) - ®2(2 PioP 12 + Pll^ + 

^‘^ 01^21 ^ ^^ 02 ^ 20 ) " 7 -~I ^0 

. ^PoiPioPu) -*■ P 

7+1 7-1 


^ 3 ^ ^(^01 ^20 ^10 *^02 


27 „ 
7 - 1 


"1 


7 + Iq 7 p + 

Poo P12 


^ 2 ' Poo 


■( 


PolPll '*" P10P02 


) '+ 3 C 


7+1 2 
/ Pnn P 


01 '^10 


T|r + 
^10 


(^/^^ POO^PU ^ 2 ^ ^ POO ^ ^ PloPoi )"*"!! 

^/■'^ P 00 '' P 10^12 ^ (""l 


7+1. 7. * ^ /n 7+1 7. ■ ^ C 7+1 7-l_ 2\^ 

Poo P02 2 Poo Poi J 20 


C -^p f + P ' 

1 Poo Poi 21 Poo 22 


27 


7-1 "2[('^l 

» Oo ’'''^ Ko ' I ' i 2 + * U^j 


7+1 7 

Poo P02 
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(^30) ^11^31 ■*■ ®1^30 “^21 ■ ^^20^30 “■ ^2^10^20 


y-\ ^o( 


^3 ^00 *^10 / 


=2’’^^'’oo''V)*10 ^ 


^1 ^00 ^10^20 *^00 ^30 


7 ^ “ 2 («/^Vp 10 C - 


Pq ^’*’^Mf llf ^ _ ^ 2 . Q^ 0 ^‘*‘^ljf 

^^00 ^0^30/ .7-1 3^00 1 


(Bqo) ^^00^12 ~ ^11^01 ~ ®^00^20 ^10^10 " 7 - 1 ^00 


(Bqi) ^^ 00^13 ^^11*^02 ” ® (*^00^21 ^ 01 ^ 20 ) ■*■ 


1 n 7+2„ 7+1. 


^(^10^11 ■*■ ^11^10) ■ 7 - 1 ^1 ^ 


00 ^01 


(B02) 12poo^ll<. - 31 ^iiPo 3 + 2pQ2%2 " 3 Poi% 3 " ®(p00^22 P01%L ■*■ 

*^ 02 ^ 20 ) ■'^ ^(^10^12 '^ll^U *^ 12 ^ 10 ) 


Poo ^02 *^2 Poo ^01 
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3^^ 


(B03) ^^00^15 " ^^11^01)- ^ ^ 03’*^12 “ ® (*^ 00^23 ^01^22 ■*■ 

^02^21 '^03^^20) ^(^lo’*^13 ‘^11^12 ^12^11 

*^13^10) ” 00 ^ ^ 01^02 


7 + 2 ^ 7-1 3 \ 

PQl ) 


r ' n ^ 
1-3 Pqo 


(Bq^) 30Poo^l6 = ^'^uPo^ ^^Ok'^12. 3Po3^13 - ^02^14 " ^5Poi^15 - 

® (‘^00^24 ^01^23 *^02^22 *^03^21 ^04^20 ) 

^(^lO^lll- ■*■ ^11^13 ^12^12 '^13^11 

y 3 1 [^I^^^POO^^^Pq!)- ■*■ ^2^^^P00^(p02^ ^^01^03) 

3C3^-^^POo’'"^P01^P02 + 



^*^00^22 


“ ^01 ^21 ■*■ ^^10^12 ^11^11 " ^^00^30 ^^^10^20 

r. ... , ^®2 7+2.,. ^^1 _ 7+l_ 

°^20^10 " ®10^0 ” 7-1 *^00 ^10 ”7-1 ^00 ^10 
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^^ 00^23 “ ^PQ2^21 ^.^11^12 ■’■ (^^12 ®P2o)^H “ 

^‘‘(‘>00*31 ^ ‘'oi»3o) "®(''l0+21 * fll*22) * ®‘>a'»10 - 

(®io^n ■*■ Sn^io) " y _ i('^oo^ ’*'11 ■*■ 


(®12) 

^^^00^24 "" "^Poi^23 (^^02 ^ ^^*^10)^22 ■*■ ( 

,^Po3 '*’ ^^*^11)^21 


^^^12’*^20 ^'^11^13 

V 

(6pi2 + 8p2o)'*^12 


(3Pi3 + 8p2i)^ll + ®P22^10 " 

[•^00^32 ^01^31 


‘’o2'*'3o) ■ (®10*h.2 * ®ll\l ■’■ ®12*\o) ■ 


7 - 1[ 

(c/-= 

7fti r 



POO^^^'^12 ■*■ ^1^^^P00^'*^^P01^11 

p 2 + c ^■*■2 7+lp U 

^00 PqI 1 ^00 ^02ri0 

Poo^^^Pi2 ■•■ ‘^i^^^Poo^PqiPii ■*■ 

7-I ■ 2 7+1 7 \ 

Poo Poi ^1 ' Poo P02/P10 


+ 
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(®2o) ^^00^32 " ^^10^22 ^2l(‘^U ~ ^^01®lo) “ ^^12 (*^00^10^ " ^2o) 

^u(^®10^11 " ^^10^^01 " ^21 ^^01®2o) ^31^01 ■ 

l^j^^l^Poo ■ ^t^oP^o - 2^20 ('^10^10 " ^20) 

^lo(^® 10^20 ” ~ ‘^ 3 o)J 

^1 ! y y+ 2 . y y 2\ 

y _ 2_{^1 PqO ^20 ^2 PqO PlO ) ~ 


^—(0 + C Tlf ] - p ^ 

- 1\'^00 ^20 1 ^00 ^ 10 ^ 10 / 7-1 ^00 ^10 


In the above expressions, a number of symbols are used for abbreviation. 
They are defined as follows: The guj^'s are coefficients of the 

quantity Pjn/Po : 



Hence 


SjnO 


^mO 

P(X) 


_ P ml PmOPOl 

’™1 Pnn p ^ 
00 


Sm2 


Pm2 

*^00 


PmlPQl ^ Pmo/PQI^ 
^00^ ^oox'^oo^ 
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The Cn™'s are the binomial coefficients sucb that 


'n 


that Is, 


m 


(1 + x)“ = ^ 
n=0 


m(m - 1 ) . . . (m - n + 1) 


n 


Finally the Djjj's and h^^'s are simply abbreviations: 


Dl = Cl ^^900^ - ^POO 


2ao^ 


7+1q 7-1 

^ _ 1 "2 Poo 


D 2 = - 20 


"^10 = 2Poi^u + 


^11 “ ^^oA 2 ■*■ (^02 ® ^10)^11 ®^u\o 


^12 ^‘^ Ol^ia (®^02 ^ ® Plo )\2 (^^03 ^*^ 12^0 


For an expemslon starting from the shock, the conditions at the 
shock are given by Identifying equations (12) and (l3b) vlth equation (20). 
The equation obtained from equating the coefficients of a certain power 
of X In the condition for p, say. Is now to be denoted by a super- 

( 2) 

script to p. Thus p' stands for the equation resulting from the 
2 

coefficients of x In the condition for p; similarly, i(r 


stands 
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for "the same in 'tlie condition for ''|f. Then the following are found to 
hold: 



^01 ^10^1 "^1^00^1 


^02 ^ 1*^11 ■*' ^ 1^^20 "^ 2*^10 ^ oof ‘^ 1 ^ 2^1 ^ 

P2I '*' ^l^pso “ "P03 " ^1^12 " ^2Pn " ^sPlO " ^^1^2P20 '*' 

Poo[}pl^3^1 ^ ■ ^2^2)^! ^ ■ 

(ci3 - - 3C3^,^3 + 8c3^o2)^^-3] 

(♦< 2 )) ’I'u + = 0 

O 

^ 21^1 ■*■ ^ 0 ^ 1 ^ = "'^12 " ^ 20^2 
(^fW) t3i3j_2 + = _3}r^^ _ - ^21^2 - ^22^1 ~ 


One may proceed in the order Indicated helow to arrive at all the 
coefficients • • •/ Pq^* The right-hand side of each 

expression is the irnknown which may he determined from the left-hand 
side in terms of quantities already obtained. 
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(^0) 


identity between and 


PQO 




simultaneously 






1(0 
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As a result, for a fourth-degree expansion, one only has to go as 
far as equation solving a total of 16 relations. For a sixth- 

degree one, there are altogether 30 relations. The increase of amount 
of work gives an indication as to what must he expected for a finer 
approx -? ma t ion . 

Next, for an expansion starting from the stagnation point, the same 
set of equations [Aqq^ to ^B2 q^ still holds. The quantities will he 

denoted by D^, and so forth. The conditions resulting from 

eqxiation (23) take the same form as . . ., except -that 

the 3's are to be replaced by the e's. But the conditions resulting 
from the densl-by -variation along the body are missing. The pre-vious 
order for the successive solution of the unkno-wns may still be followed, 
while one arbitrary parameter must now be Introduced in each set of the 

simultaneous equations like the set ^-A^q)^ and . Thus, for 

a sixth-degree expansion one may take ^13^ ftnd as the param- 

eters and break up the simultaneous equations into separate ones. Mean- 
while the fact liCj^Q = 0 greatly siBq>lifies -the coaputation. The neces- 
sary steps, in fact, are reduced to: 



(^03) ^ ^03 


(®0l) ^ ^^21 
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h2 


(All)" 

-> ppi 


— > ^30 

(■*■ 20 ) " 

^20 

(®io) - 

-^^22 

(®02)“ 

-^^11, 

(^ 05 ) ~ 

^05 




Hence a total of only 8 steps are sufficient for a fourth-degree 
expansion, and.l6 steps, for a sixth-degree one. Further, there is no 
longer the need to solve some of the equations simultaneously. 
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